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In this paper, we show how the method of field theoretical renormalization group may be used to analyze 
universal shape properties of long polymer chains in porous environment. So far such analytical calcula- 
tions were primarily focussed on the scaling exponents that govern conformational properties of polymer 
macromolecules. However, there are other observables that along with the scaling exponents are universal 
(i.e. independent of the chemical structure of macromolecules and of the solvent) and may be analyzed within 
the renormalization group approach. Here, we address the question of shape which is acquired by the long 
flexible polymer macromolecule when it is immersed in a solvent in the presence of a porous environment. 
This question is of relevance for understanding of the behavior of macromolecules in colloidal solutions, near 
microporous membranes, and in cellular environment. To this end, we consider a previously suggested model 
of polymers in d-dimensions [V. Blavats'ka, C. von Ferber, Yu. Holovatch, Phys. Rev. E, 2001, 64, 041102] in 
an environment with structural obstacles, characterized by a pair correlation function h(r), that decays with 
distance r according to a power law: h(r) ~ r~ a . We apply the field-theoretical renormalization group ap- 
proach and estimate the size ratio (Rc)/(F%) and the asphericity ratio A d up to the first order of a double 
e = 4— cf, 8 = 4— a expansion. 
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1. Introduction 

Polymer theory belongs to and uses methods of different fields of science: physics, physical 
chemistry, chemistry, and material science being the principal ones. Historically, the structure 
of polymers remained under controversial discussion until, under the influence of the work by 
Hermann Staudinger the idea of long chain- like molecules became generally accepted. In this 
paper we will concentrate on the universal properties of long polymer chains immersed in a good 
solvent, i.e. the properties that do not depend on the chemical structure of macromolecules and of 
the solvent. Usually a self-avoiding walk model is used to analyze such properties [2, 3]. At first 
glance such a model is a rough caricature of a polymer macromolecule since out of its numerous 
inherent features it takes into account only its connectivity and the excluded volume modeled by 
a delta-like self-avoidance condition. However, it is widely recognized by now that the universal 
conformational properties of polymer macromolecules are perfectly described by the model of self- 
avoiding walks. It is instructive to note that the idea to describe polymers in terms of statistical 
mechanics appeared already in early 30-ies due to Werner Kuhn [4( and already then enabled 
an understanding and a qualitative description of their properties. The present success in their 
analytic description which has lead to accurate quantitative results is to a large extent due to the 
application of field theoretic methods. In the pioneering papers by Pierre Gilles de Gennes and his 
school [2] an analogy was shown between the universal behaviour of spin systems near the critical 
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point and the behaviour of long polymer macromolecules in a good solvent. In turn, this made it 
possible to apply the methods of field theoretical renormalization group Q to polymer theory. 

In spite of its success in explaining the universal properties of polymer macromolecules, the 
self-avoiding walk model does not encompass a variety of other polymer features. Different models 
and different methods are used for these purposes. In the context of this Festschrift it is appropriate 
to mention the approach based on the integral-equation techniques which is actively developed by 
Yurij Kalyuzhnyi and his numerous colleagues [6|. In particular, this approach has enabled an 
analytic description of chemically associating fluids and the representation of the most important 
generic properties of certain classes of associating fluids It is our pleasure and honor to write 
a paper dedicated to Yurij Kalyuzhnyi on the occasion on his 60th birthday and doing so to wish 
him many more years of fruitful scientific activity and to acknowledge our numerous common 
experiences in physics and not only therein. 

In this paper we will analyze the shapes of polymer macromolecules in a good solvent. Flexible 
polymer macromolecules in dilute solutions form crumpled coils with a global shape, which greatly 
differs from spherical symmetry and is surprisingly anisotropic, as it has been found experimentally 
and confirmed in many analytical and numerical investigations (il-l22|. Topological properties of 
macromolecules, such as the shape and size of a typical polymer chain configuration, are of interest 
in various respects. The shape of proteins affects their folding dynamics and motion in a cell 
and is relevant in comprehending complex cellular phenomena, such as catalytic activity [23]. 
The hydrodynamics of polymer fluids is essentially affected by the size and shape of individual 
macromolecules [24j|; the polymer shape plays an important role in determining its molecular 
weight in gel filtration chromatography [25]. Below we will show how the shape can be quantified 
within universal characteristics and how to calculate these characteristics analytically. 

The set up of this paper is as follows. In the next section we present some details of the first an- 
alytical attempt to study the shape of linear polymers in good solvent, performed by Kuhn in 1934. 
Since then, the study of topological properties of polymer macromolecules was developed, based on 
a mathematical description, which is presented in section 3 along with a short review of the known 
results for shape characteristics of flexible polymers. In section 4, we present details of the appli- 
cation of the field-theoretical renormalization group approach to the study of universal polymer 
shape characteristics. Section 5 concerns the effect of structural disorder in the environment on the 
universal properties of polymer macromolecules. The model with long-range-correlated quenched 
defects is exploited, and the shape characteristics are estimated in a field-theoretical approach. We 
close by giving conclusions and an outlook. 



2. Shape of a flexible polymer: Kuhn's intuitive approach 

The subject of primary interest in this paper will be the shape of a long flexible polymer 
macromolecule. That is, assuming that a polymer coil constitutes of a large sequence of monomers, 
does this coil resemble a globe (which would be the naive expectation taken that each monomer 
is attached at random) or does its shape possess anisotropy. And, if it is anisotropic, what are the 
observables to describe it quantitatively? In this section we start our analysis with a very simple 
model that allows us to make some quantitative conclusions. The aim of the calculations given 
below is to explain how the anisotropy of the polymer shape arises already within a random walk 
model. That is to show that the anisotropy is essentially not an excluded volume effect (although 
its strength is effected by the excluded volume interaction as we will show in this paper) but rather 
it is an intrinsic property arising from random walk statistics. 

The analysis given below is inspired by Kuhn's seminal paper [4|. However, Kuhn's explanation 
is based on combinatorial analysis and application of Stirling formula, while here we suggest a 
derivation based on the application of Bayesian probability [26]. Let us consider the so-called 
Gaussian freely jointed chain model consisting of N connected bonds capable of pointing in any 
direction independently of each other. Any typical configuration of such a chain can be represented 
by the set of bond vectors {a n }, n = 1, . . . , N, such that: 

(S 2 )=a 2 =d£ 2 , (2.1) 
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(b) 



Figure 1. a) Schematic presentation of a polymer chain with the end-to-end distance (position 
vector of iVth monomer) denoted as r*i . The component of the position vector of iV/2th monomer 
in direction perpendicular to r\ is denoted by r*2 • Vector r*3 denotes the component of position 
vector of iV/4th monomer in direction perpendicular to both r\ and r"5 . b) In the Edwards 
continuous chain model, the polymer is represented by a path f(t) parameterized by ^ t ^ T 
(see section 4). 



here d is the space dimension and the angular brackets stand for averaging with respect to different 
possible orientations of each bond. Fixing the starting point of the chain at the origin, one gets for 
the end-to-end distance R c and its mean square: 



JV 



(iil) = Na2 



(2.2) 



To get the second relation, one has to take into account that (didj) — for i =/= j, since random 
variables Si are uncorrelated. Due to the central limit theorem, the distribution function of the 
random variable R c — \R C \ takes on a Gaussian form: 



P{Re) 



(2vr(i?2)) 



d/2 



(2.3) 



Numerical factors in (|2.3p can be readily obtained from the normalization conditions for the distri- 
bution function and its second moment. Let us consider the three-dimensional (d — 3) continuous 
chain, such that 



3 /2 r 2 

e"l^47ri?Mi?o 



(2.4) 



defines the probability that the end point of the chain is located in the interval between R c and 
R e + dR c . Following Kuhn, let us take the mean value of R e as one of the shape characteristics of 
a chain and denote it by ¥\ = (R e )- One gets for its value: 



OO . 

/ R e P(R e )dR e =2iJ—. 



(2.5) 



To introduce two more observables that will characterize the shape of the chain, let us proceed 
as follows. Having defined the end-to-end vector R e , let us point the z-axis along this vector (see 
figure [1] (a)). Now, since both the starting and the end points of the chain belong to the z-axis, the 
projection of the polymer coil on the ccy-plane has a form of a loop. It is natural to expect that 
the largest deviation of a point on this loop from the origin corresponds to the iV/2th step. Let 
us find the distance in the plane from the origin to this point (we will denote it by r-i hereafter) 
and take it as another shape characteristics of the chain. Note, that vector f*2 lies in the xy-plain 
and therefore is two-dimensional. To do so, we have to find the distribution function of a position 
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vector of a point on a loop. It is convenient to find such a distribution using the Bayes theorem 
for conditional probability [26]. The theorem relates the conditional and marginal probabilities of 
events A and B, provided that the probability of B does not equal zero: 

In our case, P(A\B) = -P2M is the probability that the coordinate of the chain after N/2th step is 
given by a (now two-dimensional) vector r under the condition, that after N steps its coordinate is 
r = 0. Then, P{B\A) is the probability for the chain that begins at the point with the coordinate 
r after N/2 steps to return to the origin. Correspondingly, the prior probability P(A) is the 
probability that the coordinate of the chain after AT/ 2th step is given by a vector r (it is the so- 
called "unconditional" or "marginal" probability of A, in our case it is given by equation (]2.3p for 
d = 2). P(B) is the prior or marginal probability of B, in our case this is the probability for the 
chain that starts at the origin to return back in N steps, i.e. a probability of a loop of N steps. 
Realizing that for our case P{A) = P(B\A), that is probabilities to reach point r starting from the 
origin is equal to the probability to reach the origin starting from the point r we get: 

P 2 (r) = P(A) 2 /P(B), (2.7) 

where P(A) is given by equation (12. 3p for d = 2 and P(B) may be found from the normalization 
condition. In the continuous chain representation we get for the probability that the AT/ 2th bond 
of the chain is located in the interval between r and r + dr: 

P2(r)dr = 2*N/2P P ( B f ~^ e ~^ 27rrdr - ( 2 - 8 ) 
The mean value follows: 

00 00 
r 2 = J r P 2 {r)dr = j^ J e~^ r 2 dr = ^^. (2.9) 


Now, let us point an x-axis along fa (again see figure [1] (a)) and repeat the above reasonings 
concerning the maximal distance in y-coordinate, i.e. concerning the y-coordinate of the AT/4th step 
of the chain. We will denote it by r$ . The result readily follows by the analogy with equation (|2.9p 
taking into account that now the radius- vector is one-dimensional: 

00 00 . 

f ( 16 \V 2 f _4r 2 i In 

r, s yrfl,(r)dr=^j J e ^ r dr = -^-. (2.10) 



Comparing (12. 5p , (|2.9p , and (|2.10p one concludes that the shape of the chain is characterized 
by three unequal sizes r% , F2 , and ¥3 with the following relations: 

^i = -~2.55, 'Q- = 4V2 ~ 5.66. (2.11) 

The above relations were first obtained by Kuhn [J] and lead to the conclusion that a polymer chain 
even if considered as a chain of mutually intersecting steps (i.e. without account of an excluded 
volume effect) does not have a shape of a sphere but rather resembles an ellipsoid with unequal 
axetQ. To check how this prediction holds, we have performed numerical simulations of random 
walks on simple cubic lattices, constructing trajectories with the number of steps TV up to 400 
and performing the averaging over 10 6 configurations. As one can see from figure [5] the results of 
simulations are in perfect agreement with the data of (|2.1ip . 




1 As it was stated in Kuhn's paper, the most probable shape of a polymer is a bend ellipsoid, of a bean-like shape: 
. . verbogenes Ellipsoids (etwa die Form einer Bohne). . ." 0|. 
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Figure 2. Ratios ri/ri (a) and n/rs (b) as functions of the chain length, results of computer 
simulations. Analytic estimates give: T\jr% — 2.55 and f\jfz ~ 5.66 (see equations (|2.5p . (|2.9p , 
(l2~T0l)l. 



Here it is worth mentioning another far going prediction of the same paper [4|. Discussing 
how might the excluded volume effect the polymer size, Kuhn arrives at the relation between 
the mean square end-to-end distance and the number of monomers which in our notations reads: 

\j (R%) = iN v with v = 1/2 + e. Although this result is suggested on purely phenomenological 
grounds, its amazing feature is that the power-law form of the dependence is correctly predicted 
(cf. equation (|3.4p from the forthcoming section). Moreover, Kuhn has estimated e by considering 
the excluded volume effect for a 5-segment chain for which he found an increase of 20%, a result 
we have verified by exact enumerations. Assuming the power law form, e is estimated as: e « 0.11 
and thus v w 0.61, which is perfectly confirmed later, e.g. by Flory theory 0|, which in d = 3 gives 
v = 3/5. The notation for the correction used by Kuhn by coincidence is the same as that used 
much later in the famous e-expansion [H| to develop a perturbation theory for calculation of this 
power law by means of the renormalization group technique (see equation (|4.12p for the first order 
result). Before explaining how this theory is applied to calculate polymer shape characteristics, let 
us introduce observables in terms of which such description is performed. 



3. Description of polymer shape in terms of gyration tensor and combina- 
tions of its components 

Let R n — {x 1 ^ ,...,x„} be the position vector of the nth monomer of a polymer chain (n = 
1, . . . , N). The mean square of the end-to-end distance R c of a chain thus reads: 

(i? c 2 ) = (|i? N -i? 1 | 2 ), (3.1) 

here and below, (. . .) denotes the averaging over the ensemble of all possible polymer chain con- 
figurations. The basic shape properties of a specified spatial conformation of the chain can be 
characterized [8|, Q in terms of the gyration tensor Q with components: 

1 N 

Qij = — x CM)( X n ~ X Cm)> ^1 3 = 1, • • • , <i, (3-2) 

71=1 

with Xq M = Yln=i X n/N being the coordinates of the center-of-mass position vector Rqm ■ 

The spread in the eigenvalues Ai of the gyration tensor describes the distribution of monomers 
inside the polymer coil and thus measures the asymmetry of the molecule; in particular, for a 
symmetric (spherical) configuration all the eigenvalues Ai are equal, whereas for the so-called 
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Figure 3. Seen from far away, polymer coil may resemble the objects of different from. Here, we 
distinguish sphere-like (a), prolate (b), and oblate (c) conformations. Eigenvalues of correspond- 
ing gyration tensor (|3.2|) satisfy: Ai ~ A2 ~ A3 (a), Ai > A2 ~ A3 (b), and Ai ~ A2 > A3 (c). 
Correspondingly, asphericity (|3.6ll and prolateness 1)3. 7[) of these conformations satisfy: Ad = 0, 
5 = (a), < Ad < 1, < S < 2 (b), < A d < 1, -1/4 < S < (c). See the text for more 
details. 



prolate and oblate configurations in d — 3 (see figure [3]) the eigenvalues satisfy Ai 3> A2 ~ A3 
and Ai w A2 > A3 correspondingly. Sole and Stockmayer [8| introduced the normalized average 
eigenvalues Ai of the gyration tensor as a shape measure of macromolecules. Numerical simulations 
in d = 3 dimensions give {(Ai), (A2), (A3)}={0. 790, 0.161, 0.054} [16j, indicating a high anisotropy 
of typical polymer configurations compared with the purely isotropic case {1/3, 1/3, 1/3}. 

While in simulations the eigenvalues of the gyration tensor Q can easily be calculated and 
averaged, different invariants have been devised for theoretical calculations. As far as Q has three 
eigenvalues in d = 3, one may construct three independent combinations of invariants. These are 
the square radius of gyration R G , the asphericity Ad and the prolateness S as elaborated in the 
following. The first invariant of Q is the squared radius of gyration 

N d 

i?G = ^E( i? »- i? CM) 2 = J]g«-TrQ, (3.3) 

n— 1 i— 1 

which measures the distribution of monomers with respect to the center of mass. To characterize 
the size measure of a single flexible polymer chain, one usually considers the mean-squared end- 
to-end distance (R 2 ) (|3.ip and radius of gyration (R%) (13. 3p . both governed by the same scaling 
law: 

(R 2 ) ~ (R 2 G ) ~ N 2 \ (3.4) 

where N is the mass of the macromolecule (number of monomers in a polymer chain) and v is 
a universal exponent (v > 1/2 (d < 4), v = 1/2 (d 4)). The ratio of these two characteristic 
distances, the so-called size ratio: 

g d =(Rl)/(R 2 G ), (3.5) 

appears to be a universal, rotationally-invariant quantity (gd > 6 (d < 4), gd = 6 (d ^ 4)) (l5| . 

Let A = Tr Q/d be the mean eigenvalue of the gyration tensor. Then one may characterize the 
extent of asphericity of a polymer chain configuration by the quantity Ad defined as [Io| : 

^-did-l)^ t "<*-l(TrQ) 2 ' ^ 

with Q = Q — X I (here I is the unity matrix) . This universal quantity equals zero for a spherical 
configuration, where all the eigenvalues are equal, A; = A, and takes a maximum value of one in 
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the case of a rod-like configuration, where all the eigenvalues equal zero except one. Thus, the 
inequality holds: ^ Ad $C 1 . Another rotational invariant quantity, defined in three dimensions, 
is the so-called prolateness S: 



nti(^-A) detQ 
A 3 (TrQ) 3 



5= ^=1^-^ =27^^. (3.7) 



If the polymer is absolutely prolate, rod-like (Ai ^ 0, A2 = A3 =0), it is easy to see that S 
equals two. For absolutely oblate, disk- like conformations (Ai = A2 , A3 = 0), this quantity takes 
on a value of —1/4, while for a spherical configuration 5* = 0. In general, S is positive for prolate 
ellipsoid- like polymer conformations (Ai 3> A2 ~ A3) and negative for oblate ones (Ai ~ A2 ^> A3), 
whereas its magnitude measures how much oblate or prolate the polymer is. Note that since A and 
the quantities in ()3.3[) (|3.7p are expressed in terms of rotational invariants, there is no need to 
explicitly determine the eigenvalues Ai which greatly simplifies the calculations. 

The average of the quantities (|3.3|) - (|3.7p for a given polymer chain length N, denoted as (...), 
is performed over an ensemble of possible configurations of a chain. Note that some analytical 
and numerical approaches avoid the averaging of the ratio in (|3.6|) , (|3.7p and evaluate the ratio of 
averages: 

Ad = 1 ±^i, i.nLM, ( 3 . 8 ) 

<<(<<- i)tr (A-> (A s > 

which should be distinguished from the averaged asphericity and prolateness: 

Contrary to (Ad) and (S), the quantities (|3.8p have no direct relation to the probability distribution 
of the shape parameters Ad and S. As pointed out by Cannon [11], this definition overestimates the 
effect of larger polymer configurations on the mean shape properties and suppresses the effect of 
compact ones. This artificially leads to overestimated values for shape parameters. The difference 
between (Ad) and Ad was found to be great (see table [IJ. 




Table 1 . Size ratio, averaged asphericity and prolateness of flexible polymer chains on regular 
two- and three-dimensional lattices. MC: Monte Carlo simulations, DR: direct renormalization 
approaches, a: reference [l3| . b: (T3 |. c: [l7| . d: [ljj, e: 



d 


Method 


9d 


(A d ) 


A d 


(S) 


S 


2 


MC 


7.14±0.03 a 


0.501 ± 0.003 6 


0.625 ±0.008 b 






2 


DR 


7.003 c 




0.558 d 




0.899 d 


3 


MC 


6.249 ±0.03 e 


0.431 ± 0.002 e 


0.546 ±0.008 6 


0.541 ± 0.004 e 




3 


DR 


6.258 c 


0.415 6 


0.529 d 




0.893 d 



Numerous experimental studies indicate that a typical flexible polymer chain in good solvent 
takes on the shape of an elongated, prolate ellipsoid, similar to what was shown in section [3] for a 
random walk. In particular, using the data of x-ray crystallography and cryo-electron microscopy, 
it was found that the majority of no n-g lobular proteins are characterized by A3 values from 0.5 



to 0.7 and S values from to 0.6 [21] , |22j. The shape parameters of polymers were analyzed 
analytically, based on the direct renormalization group approach fiol . Hil . QjJ, and estimated in 
numerical simulations 0, [13, EH llij l- Table [T] gives typical data for the above introduced shape 
characteristics of long flexible polymer chains in d = 2 and d = 3. 
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Since the above shape and size characteristics of polymer macromolecules are universal, i.e. 
independent of the details of their chemical structure, they were (along with polymer scaling 
exponents) the subject of analysis by field-theoretical renormalization group approaches. In the 
subsequent section we will introduce this approach as used to calculate polymer shapes. 



4. Field-theoretical renormalization group approach to define 
polymer shape 

Aronovitz and Nelson [ic] | developed a scheme, allowing one to compute the universal shape 
parameters of long flexible polymers within the frames of advanced field theory methods. 

Here, we start with the Edwards continuous chain model [27| . representing the polymer chain 
by a path f(t), parameterized by ^ t ^ T (see figur^U (b)). The system can be described by the 
effective Hamiltonian H: 



The first term in (|4.1j) represents the chain connectivity, whereas the second term describes the 
short range excluded volume interaction with coupling constant Uq . 

In this scheme, the gyration tensor components (equation (|3.2|Q can be rewritten as: 



T 

1 

2T2 



dti / dt 2 [r^O-r^)] [^(tij-r'te)], i,j = l,...,d. (4.2) 



o 



The model (I4.1j) may be mapped to a field theory by a Laplace transform from the Gaussian 
surface T to the conjugated chemical potential variable (mass) ^iq according to [28, 29]: 

i( M o) = J dTcxp[- Mo T]Z(T), (4.3) 

where Z(T) = f D[r] exp(— T-L) is the partition function of the system as function of the Gaussian 
surface and J D[r] means an integration over all possible path configurations [30]. Exploiting the 
analogy between the polymer problem and 0(m) symmetric field theory in the limit m — > 
(de Gennes limit) it was shown [28| that the partition function of the polymer system is related 
to the m — 0-component field theory with an effective Lagrangean: 



1 l 217"/ m2 i |Y7l7 \|2\ 1 U (Tit \\2 



(4.4) 



C = J d d x 

Here, (j> is an m-component vector field <f> = {tjr, . . . , 4> m ) and: 

Z( Po ) = f D&]e- C . (4.5) 



One of the ways of extracting the scaling behavior of the model (I4.4|) . is to apply the field- 
theoretical renormalization group (RG) method @] in the massive scheme, with the Green's func- 
tions renormalized at non-zero mass and zero external momenta. The Green's function Gq^'^ can 
be defined as an average of N field components and L ^-insertions performed with the corre- 
sponding effective Lagrangean C: 

§ CE h+ ]>>) G o W ' L) ({ fc }; M; m§ ; «o) = 

3 i(fc 1 H 1+W rj}^ (ri j _ _ _ ^( ri )^ (i?1 ) . . . <j>(R N ))£°*d d Ri . . . d d R N d d r x . . . d d r L , (4.6) 
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here {k} = (k\ , . . . , ki,), {p} — (pi , . . . ,pn) are the sets of external momenta and Ao is a cut- 
off [3l|. The renormalized Green's functions G^' L ' are expressed in terms of the bare vertex 
functions as follows: 

G^' i) ({*};{p};/&«o) = z^ /2 z~ 2 L G^' L \{kY{p}-^u), (4.7) 

where Z^ , Z^i are the renormalizing factors, /i, u are the renormalized mass and couplings. 

The change of coupling constant uq under renormalization defines a flow in parametric space, 
governed by the corresponding /3-function: 

where I is the rescaling factor and | Q stands for evaluation at fixed bare parameters. The fixed 
points (FP) of the RG transformation are given by the zeroes of the f3- function. The stable FP u*, 
corresponding to the critical point of the system, is defined as the fixed point where -^j^ \u=u* 
has a positive real part. The flow of renormalizing factors , Z^i defines the RG functions 7</,(m), 
7^2 (it). These functions, evaluated at the stable accessible FP, allow us to estimate the critical 
exponents. 

Exploiting the perturbation theory expansion in parameter e = 4 — d (deviation of space 
dimension from the upper critical one), one receives within the above described scheme up to the 
first order in e the well-known results for the fixed points: 

w rw = ■ stable for e < 0, (4.9) 

3e 

"saw — ~4~ ! stable for e > 0. (4-10) 

Here, urw describes the case of simple random walks (idealized polymer chain without any in- 
trachain interactions), and wsaw is the fixed point, governing the scaling of self-avoiding random 
walks. 

Evaluating the RG functions ^y^iu) and 7^2 (u) at the above fixed points, one gets the familiar 
first-order results (see e.g. Q) for the critical exponents v and 7, that govern the scaling of the 
polymer mean size (|3.4[) and the number of configurations, correspondingly: 

"hw = g > 7kw = 1 , (4.11) 

1 £ £ 

z^saw = 2 + Jq ' 7sAW = 1+ 8 ' ( 4 ' 12 ' ) 

Following reference the averaged moments of gyration tensor Q (|4.2I) . which are needed 
to determine the polymer shape characteristics (|3.3[) . (|3.6p . (|3.7p . can be expressed in terms of 
renormalized connected Green's functions (|4.7I) . in particular: 



2\2XJ r(7+2w-2)Gg)(o,o,u') 



Here, the following notations are used: 

r =11 

13 dq i dqi 



q=Q 

d d d d 
im dq\ dq[ dq% dq\ 



2Xj r( 7 +4 J /+4) G ( i 2) (0j0iU * ) 
G^ 2) (0,0; q,-q;u*), (4.15) 



G < i A \0 1 0;q 1 ,~q 1 ,q 2 ^q 2 ;u*), (4.16) 

<;=0 



X is non-universal quantity, v and 7 are the critical exponents, d/dq\ means differentiation by the 
i-component of vector q\ , G|? and G|f are the renormalized connected Green's function with 
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(a) W (c) 



(d) (e) 

Figure 4. Contributions to the Green's function G' 2 ' 2 ' up to one-loop level. Solid lines denote 
propagators fi 2 +k 2 , wavy lines illustrate the insertions of the type (p 2 , loops imply integration 
over internal momenta. 

2 external legs and 2 insertions ip 2 (q), (p 2 (—q) and 4 insertions ip 2 (qi), ip 2 (~qi), <^ 2 (92) , ( f 2 (~q2) 
respectively, calculated at the fixed point for zero external momenta. 

The isotropy of the original theory implies, in particular, that (TrQ) = d(Qu), so that: 

(R 2 G )=d(Q u ). (4.17) 

The mean-squared end-to-end distance (R 2 ) can be expressed as: 

r( 7 ) (v 2 fe 4 2) (fc,-fc,u*)] 



fe=0 



2XJ T( 1+ 2v) 4 2) (0,0,u*) 



(4.18) 



where V 2 means differentiation over components of external momentum k. 

One can easily convince oneself, that not-universal quantities cancel when the ratio (|3.5[) is 
considered: 

' 72^7(2)/ 



_ (R 2 ) _ T( 7 +2^+2) (V fc GV( fc .-M* 

9 — IT>2 • 



fc=0 



( R %) r (7+2^) (v?Gg" 2) (0,0; g ,-g;«*)) 



(4.19) 



g=0 



In derivation of (|4.19p we made use of an obvious relation J2i=i Gu — V 2 G^ 2 ' 2 ' ) . 

Computing the asphericity we follow equation (|3.8p considering Ad as the ratio of averages: 

a d ( Tr Q 2 ) um 

Ad - d-l((TrQ) 2 )' (4 ' 20) 

This definition allows one to directly apply the renormalization group scheme described above, and 
express the Ad in terms of the averaged moments of the gyration tensor (|4.13|) , (|4.14|) : 

_ (Q^+djQp - (QuQjj) 
Ad ~ (Q|)+d( d -l)(Q i ,Q,,) > ^ 

One can again easily convince oneself that all the non- universal quantities in equations (|4. 13[) , (|4.14D 
cancel when calculating (|4.21[) . Note that within the RG approach we will focuss only on two 
universal shape characteristics, namely the size ratio (|4.19p and asphericity (|4.2ip . as far as the 
calculation of prolateness S is particularly cumbersome. 

To estimate (|4.19p . one needs to calculate the Green's functions G^ 2,2 \ presented diagrammati- 
cally in figure [U Applying the renormalization procedure, as described above, one finds to the first 
order of the renormalized coupling u (i.e., in the so-called one-loop approximation): 

G ( i' 2) (0,0; q,-q;u) = - IJ—uhfaq) ~ \uh(Q, 0, q )+i-±-uh(0, 0), (4.22) 
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(a) 



(b) 



(c) 




M 



If) 




(h) 



(i) 



0) 



Figure 5. Contributions to the Green's function G^ 2 ' 4 - 1 up to one-loop level. Solid lines denote 
propagators fi 2 +k 2 , wavy lines illustrate the insertions of the type (p 2 , loops imply integration 
over internal momenta. 



the one-loop integrals Ii are listed in appendix A. Performing an e = 4 — ri-expansion of loop 
integrals (see appendix A for details) and differentiating over components of vector q, we found: 



(v 2 q G ( £' 2) (0,0; q,-q;u)) 



q=0 

= -2. 



19 
— i 
18 



fc=0 



(4.23) 
(4.24) 



Substituting the values of fixed points (I4.9[) and (I4.10[) into the ratio of these functions, one receives 
the corresponding value of the g-ratio: 



ffRW — 6 , 
3SAW = 6 + 



16 



(4.25) 
(4.26) 



Equation (|4.26p presents a first order correction to the (/-ratio caused by excluded volume interac- 
tions. 

To compute the averaged asphericity ratio using (|4.21l) . one needs the Green function G^ 2,4 ) 
with four insertions ip 2 (qi)/2, (p 2 (—qi)/2, (p 2 (q2)/2, ip 2 (—q2)/2, which is schematically presented in 
figure [5] Applying the renormalization scheme described above, one receives an analytic expression 
for the renormalized function G^' 4 \ given in the appendix B, equation (|B1[) . Taking derivatives 
over components of inserted vectors q\ , qi and performing e-expansions of the resulting expressions 
(see appendix A for details) we arrive at the following expansions for the functions (|4.15[) . (|4.16|) : 



c 

^ XX 


= 576-1 


4028u 
15 


(*xx\yy 


= 320H 


436u 
3 ' 


(*xy\xy 


= 128-1 


308u 
5 ' 



(4.27) 



Evaluating these relations at the fixed points (|4.9[) . (|4.10p and substituting into (I4.13[) . (|4.14p . 
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(|4.2ip one receives: 




SAW 
d 



1 

2 ' 

1 15 



(4.28) 



(4.29) 




We note, that result (|4.28[l means that even within the idealized model of simple random walks, the 
shape of a polymer chain is highly anisotropic (cf. Kuhn's picture, described in section 2). Taking 
into account the excluded volume effect makes the polymer chain more extended and aspherical: 
indeed, in three dimensions (e = 1) the above obtained quantity reads: Jsaw ^ q.53. 



5. Polymer in porous environment: Model with long-range correlated disor- 
der 



In real physical processes, one is often interested how structural obstacles (impurities) in the 
environment alter the behavior of a system. The density fluctuations of obstacles lead to a large 
spatial inhomogeneity and create pore spaces, which are often of fractal structure [32| . In polymer 
physics, it is of great importance to understand the behavior of macromolecules in the presence 
of structural disorder, e.g., in colloidal solutions [33] or near the microporous membranes (34[. In 
particular, a related problem concerns the protein folding dynamics in the cellular environment, 
which can be considered as a highly disordered environment due to the presence of a large amount 
of biochemical species, occupying up to 40% of the total volume [35j|. Structural obstacles strongly 
affect the protein folding [36j. Recently, it was realized experimentally [37fl that macromolecular 
crowding has a dramatic effect on the shape properties of proteins. 

In the language of lattice models, a disordered environment with structural obstacles can be 
considered as a lattice, where some amount of randomly chosen sites contain defects which are 
to be avoided by the polymer chain. Of particular interest is the case when the concentration 
of lattice sites allowed for the SAWs equals the critical concentration and the lattice is at the 
percolation threshold. In this regime, SAWs belong to a new universality class, the scaling law l|3.4p 
holds with exponent v Pc > vsaw [38(| . The universal shape characteristics of flexible polymers in 
disordered environments modeled by a percolating lattice were studied recently in [39]. Another 
interesting situation arises when the structural obstacles of environment display correlations on 
a mesoscopic scale (4£j- One can describe such a medium by a model of long-range-correlated 
(extended) quenched defects. This model was proposed in reference [4l| in the context of magnetic 
phase transitions. It considers defects, characterized by a pair correlation function h{r), that decays 
with a distance r according to a power law: 



at large r. This type of disorder has a direct interpretation for integer values of a; namely, the 
case a = d corresponds to point-like defects, while a = <i — 1 (a = d — 2) describes straight lines 
(planes) of impurities of random orientation. Non-integer values of a are interpreted in terms of 
impurities organized in fractal structures (42| . The effect of this type of disorder on the magnetic 
phase transitions has been a subject to numerous studies [43] . 

The impact of long-range-correlated disorder on the scaling of single polymer chains was ana- 
lyzed in our previous works [3] by means of field-theoretical renormalization group approach. In 
particular, it was shown that the correlated obstacles in environment lead to a new universality 
class with values of the polymer scaling exponents that depend on the strength of correlation ex- 
pressed by parameters a. The question about how the characteristics of shape of a flexible chain 
are effected by the presence of such a porous medium was briefly discussed by us in reference (45| . 
The details of these calculations will be presented here. 

We introduce disorder into the model (|4.4[) . by redefining $ — > (iQ+SfiQ(x), where the local 



h(r) 



r 



— a 



(5.1) 
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fluctuations 5[1q(x) obey: 



((5fi (x))) = 0, 

({S^ (x)Sfi (y)))=h(\x-y\). (5.2) 



Here, ((. . .)) denotes the average over spatially homogeneous and isotropic quenched disorder. The 
form of the pair correlation function h(r) is chosen to decay with distance according to (|5.1D . 

In order to average the free energy over different configurations of the quenched disorder we 
apply the replica method to construct an effective Lagrangean 



^dis — 



lib I A d x \^i\$ a { X )\"Mv$ a { X r)+^{$i{ X )) 

n p 

+ ]T d d xd d yh(\x-y\)^ a (x)^(y). (5.3) 

a,/3=l J 

Here, the term describing replicas coupling contains the correlation function h(r) (|5.1L Greek 
indices denote replicas and both the replica (n — » 0) and polymer (m — > 0) limits are implied. For 
small k, the Fourier-transform h(k) of (|5.ip reads: 

h(k) ~ v +w \k\ a - d . (5.4) 

Taking this into account, rewriting equation (|5.3p in momentum space variables, and recalling 
the special symmetry properties of (|5.3I) that appear for to, n — » (44J, a theory with two bare 
couplings uq , wo results. Note that for a ^ d the u>n-term is irrelevant in the RG sense and one 
restores the pure case (absence of structural disorder) . As it will be shown below, this term modifies 
the critical behaviour at a < d. We will refer to this type of disorder as long-range-correlated and 
denote by LR hereafter. 

To extract the scaling behavior of the model (I5.3[) . one applies the field-theoretical renormaliza- 
tion group method following the scheme described in a previous section, with modifications caused 
by the presence of the second coupling constant wq ■ In particular, the change of couplings uq , wq 
under renormalization defines a flow in parametric space, governed by corresponding /3-functions 
(c.f. equation (|4.8[1 ): 



du 



di 



P u {u,w) = -^-— , P w (u,w) = ^— . (5.5) 
am/ o a In/ o 

The fixed points (FPs) of the RG transformation are given by common zeroes of the /3-functions. 
The stable FP (u* , w* ) that corresponds to the critical point of the system, is defined as the fixed 
point where the stability matrix = dj3\Jd\j , i,j — 1,2 possesses eigenvalues with positive 
real parts (here, Ai = u, A2 = w). 

In our previous work [44| we have found the FP coordinates for polymers in LR disorder, which 
up to the first order ofe = 4 — d, <5 = 4 — a-expansion read: 

u rw = 0' w rw = stable for 6 < 0, S < 0, (5.6) 

3e 

"saw = T ' w saw = stable for S < e/2, (5.7) 



4 

3r5 2 „ _ 3S(e - 25) 



w lr =_ ^7 FT stable for e/2 < 5 < e. (5.8) 



2(e - 6) ' Ln 2(e - 6) 



The RW and SAW fixed points restore the corresponding cases of a polymer in pure solvent 
(cf. (|4.9p . (|4.10p ). whereas the LR fixed point reflects the effect of correlated obstacles, which 
appears to be non-trivial in certain regions of the d, a plane and to govern new scaling behaviour 
in this region. For critical exponents z^lr, 7lr governing the scaling of polymer chains in the region 
of a, d, where the effect of LR disorder is nontrivial, one obtains [44|: 

*/lr - 1/2+6/8, 7LR - l+S/4. (5.9) 
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As it was explained in the previous section, to estimate the universal size ratio (|4.19p for the 
case of polymers in long-range-correlated disorder, we calculate the Green function G^ 2 ^{u, to), 
presented diagrammatically in figure @] Now, every interactive diagram appears twice, once with 
each of the two couplings u and to, respectively. For the renormalized function we obtain up to the 
one-loop approximation: 

G^ 2) (0, 0; q, -q; u, w) = - \ [uh (0, q) - wJ x (0, q)] 

q z +l 3 q z +l 

-\ [ul 2 (0, 0, q) - w.J 2 (0, 0, q)] +\^~ [uh(0, 0) - wJt(0, 0)] . (5.10) 
3 3 q z +l 

The one-loop integrals Ii are given in appendix A. Note that Ji differs from corresponding Ii 
only by an additional factor |p| a ~ d in the numerator. Differentiating over components of vector q, 
evaluating at q — and performing double e, ^-expansions of the resulting expression, we have: 



; V^ 2) (0, 0; q, -q- u, wj) | 3=q = -4 - - (5.11) 



VtG£{k,-k)) fc=Q = -2. (5.12) 

Recalling the value of the LR fixed point (|5.8[) . and inserting this into the ratio (|4.19p of these 
functions, one arrives at an expansion for the size ratio of a SAW in the presence of long-range 
correlated disorder: 

Slr = 6+- . (5.13) 

This should be compared with the corresponding value in the pure case (|4.26l) . Let us qualitatively 
estimate the change in the size ratio g, caused by the presence of structural obstacles in three 
dimensions. Substituting directly e — 1 into (|4.26p . we have for the polymer chain in a pure 
solvent: g puTe ~ 6.06. Let us recall that the effect of long-range-correlated disorder is relevant to 
a ^ d (5 ^ e) (see e.g. explanation after equation (|5.3[) L Estimates of <7lr can be evaluated by 
direct substitution of the continuously variable parameter S into equation (|5.13p . One concludes 
that increasing the parameter S (which corresponds to an increase of disorder strength) leads to a 
corresponding increase of the <?-ratio. 

To compute the averaged asphericity in correspondence with (14.2ip , we need the Green function 
G^ 2 ' 4 ' (u,w) with four insertions ip 2 (qi)/2, ip 2 (—qi)/2, ip 2 (q2)/2. ip 2 (—q2)/2, shown diagrammati- 
cally in figure [SJ The corresponding analytic expression for the renormalized function is given in 
appendix B, equation (|L32|) . Taking derivatives of this expression with respect to the components 
of the inserted vectors q\, qi we find: 

Gxx = 576H — — (u-w), 
15 

G xx \ yy = 320+-^(u - to), 
308 

G xv]xy = 128+— (u - w). (5.14) 
At the LR fixed point (JSTSJ), from P~2"Tj) we finally have: 

Again, let us qualitatively estimate the change in Ad caused by the presence of structural obstacles 
in three dimensions. Substituting directly e = 1 into (|4.29p . we have for the pure case: A p ™ c ~ 0.53. 
Estimates of A^ R can be obtained by direct substitution of the continuously changing parameter 
S into equation (|5 . 1 5|) . An increasing strength of disorder correlations results in an increase of the 
asphericity ratio of polymers in disorder. This phenomenon is intuitively understandable if one re- 
calls the impact of the long-range-correlated disorder on the mean end-to-end distance exponent v. 
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Indeed, it has been shown in [4J| that such disorder leads to an increase of v, and subsequently, 
to the swelling of a polymer chain. Extended obstacles do not favour return trajectories and as a 
result the polymer chain becomes more elongated. In turn, such elongation leads to an increase of 
the asphericity ratio as predicted by equation (|5.15|) . 

6. Conclusions 

The universal characteristics of the average shape of polymer coil configurations in a porous 
(crowded) environment with structural obstacles have been analyzed considering the special case 
when the defects are correlated at large distances r according to a power law: h(r) ~ r~ a . Applying 
the field-theoretical RG approach, we estimate the size ratio g = (Re)/ (Rq) and averaged aspheric- 
ity ratio Ad up to the first order of a double e = 4 — d, 5 ~ A — a expansion. We have revealed that 
the presence of long-range-correlated disorder leads to an increase of both g and Ad as compared 
to their values for a polymer chain in a pure solution. Moreover, the value of the asphericity ratio 
Ad was found to be closer to the maximal value of one in presence of correlated obstacles. Thus, 
we conclude that the presence of structural obstacles in an environment enforces the polymer coil 
configurations to be less spherical. We believe that the obtained first order results indicate the 
appropriate qualitative changes in polymer shape caused by long-range-correlated environment. 
However, to get more accurate quantitative results, higher orders of perturbation theory may be 
needed. This is subject of further investigations. 
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Appendix A 

Here, we present the expressions for the loop integrals, as they appear in the Green functions 
G( 2 ' 2 ) and G^ 2 ' 4 ). We make the couplings dimensionless by redefining u — u fi d ~ 4 and w — w (i a ~ 4 ; 
therefore, the loop integrals do not explicitly contain the mass [31] : 



The loop integrals Ji (equations (15.10)) and (|B2[0 differ from corresponding Ii only by an additional 
factor |p| a ~ d in the numerators. The correspondence of the integrals to the diagrams in figure [5] is: 



(e),(/): integrals h (0, qt), h (0, q 2 ), Ji(0, g>i), Ji(0, q 2 ), 

(g) : l 2 (0, qi ,q 2 ), l 2 (0,0,qi), I 2 (0,0,<? 2 ), J 2 (0, qi ,q 2 ), J 2 (0,0,qi), J 2 (0,0,q 2 ), 

J 2 (0, qi ,q 2 ), J 2 (0,0, gi ), J 2 (0, 0,^), 

(h) : l 2 (0,0,q x ), / 2 (0,0,g 2 ), / 2 (0, qi , qi+q 2 ), W,q 2 ,qi+q 2 ), J 2 (0,0,?i), 

•/ 2 (0,0,q 2 ), J 2 (0,gi ,qi+q 2 ), J 2 (0, q 2 , qi+q 2 ), 




I 3 (k l7 k 2 ,k 3 ,k 4 ) = j ■ 




dp 



[i>+fci) 2 +i] i(p+fc 2 ) 2 +i] [(p+hr+i] [b+fc 4 ) 2 +i] ' 

f dp 



[(p+k^+i] [{p+k 2 f+\] [( P +fc 3 ) 2 +i] [{p+kiY+i] [(p+fc 5 ) 2 +i] ' 
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(i) : J 3 (0,0,gi ,q 2 ), 4(0, q\ , q 2 , 91+92), 4(0, 91 ,9i ,91+92), 4(0, 92 , 92 , 91+92), 
^3(0,0,91,92), 4(0, 91 ,92 ,91+92), 4(0, 91 ,9i ,91+92), J 3 (0, 92 ,92 ,91+92), 

(j) : 14(0,0,0,91,92), 4(0,0, 91 ,92 ,91+92), 7 4 (0, 0,91 ,91 ,91+92), 4(0, 0, 92 , 92 , 91+92), 
4(0, 0, 0, 91 , 92), J 4 (0, 0, 91 , 92 , 91+92), 4(0, 0, 91 , 9i , 91+92), 4(0, 0, 92 , 92 , 91+92)- 

In our calculations, we use the following formula to fold many denominators into one (see e.g. 
book of D. Amit in reference [5|): 

1 r(«i +...+«,, 



■ x 



..an" r(a-i) . . .r(a n ) 

gi" ■■■x„_ 1 (l-a;i-...-a: n -i) a " 
x / dan.../ da;„_i— — — r7T — , n , +mmm+an , (Al) 



[xiai+ . . . +a; n _ia n -i+(l-an- • • • -i n -i)«n] ai 



where the Feynmann variables Xi extend over the domain on + • • • + x„_i ^ 1- 
To compute the (i-dimensional integrals we apply : 



dp 



i r(d/2)r(q-d/2) a _ fc2)d/2 -« 



(A2) 

(p 2 +2kp+m 2 ) a 2 r(a) 

here dp = d d p J4/(27r) rf , where the geometrical angular factor S4 = l/(2 d-1 7r <i ' 2 r(<i/2)) is sepa- 
rated out and absorbed by redefining the coupling constant. 
As an example we present the calculation of the integral: 

(p 2 + l)[(p+9i) 2 + l][(p+9 2 ) 2 + l] 
First, we make use of formula (|A1[) to rewrite: 

1 T(3) J* dxi dx 2 



(p 2 +l)[(p+9i) 2 +l][(p+9 2 ) 2 +l] [p 2 +2p(q 1 x 1 +q 2 x 2 ) + l+x 1 q 2 +x 2 q 2 ] 3 ' 

Now one can perform the integration over p, passing to the d-dimensional polar coordinates and 
making use of the formula (|A2[) : 

dp 

[p 2 +2p(9iXi+9 2 x 2 )+l+:Ei9 2 +:E 2 9 2 J 3 

^[l+xi9?(l - x 1 )+x 2 qi(l - x 2 ) - 2 Xl x 2 qiq2] d/2 - 3 . (A4) 



2r(3) 

As a result, we are left with: 



W..,.0-ir(f)r(s-f 



1 

X 



J dxi J dx 2 [l+xi9 2 (l - i-i)+x 2 92(l - x 2 ) - 2xix 2 9i9 2 ] d/2 3 . (A5) 


To find the contributions of this integral to G xx \ xx and G xx \ yy (according to (|4. 16j) ^ . we first 
differentiate the integrand in (|A5[) over the components of vectors 91 , q 2 : 

d d d d r 2 2 id/2-3 

4x|xx = 77777777 l+an9i(l - Xi)+X29 2 (l - x 2 J - 2zi2;29i92 

dgf dgf dq% dq% q 1= q 2 =o 

= 4 (5 - 3 ) [2x 2 x 2 +x 2 a; 1 (l - z 2 )(l - xi)] - 4 ^ - 3^ [xix 2 (l - a?i)(l - x 2 )+4x 2 x 2 ] , 

d d d d r 2 2 /, \ id/2-3, 

I xx \yy = — — -^-y-^-y[l+ X iq 1 (l-x 1 )+x 2 q 2 {l- X2 )-2 Xl x 2 q 1 q 2 \ |fc=&=o 

= ( 77 - 3 ) x 2 an(l - ar 2 )(l - an) - 4 ( - - 3 ) a: 2 2n(l - an)(l - x 2 ). 
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Finally, the contributions of 7 2 (0,gi ,(72) to G xx \ xx and G xx \ yy are found by performing integration 
of I xx \ xx and I xx \ yy over x\ , X2 in (| A5[) . Results can further be evaluated either fixing the value 
of space dimension d or performing an expansion in parameter e = 4 — d. Working within the last 
approach, we obtain up to the first order of the e-expansion: 





1 1 

1„ fd\„ / d\ r . f . , 1 



2\2/V 2 J J J ,uu 4 





o r I o ) r I 3 - o ) / ^1 / da; 2 4x| ra =s (A6) 



Appendix B 

In this appendix we give the expressions for renormalized Green functions G R 2 ' 4 ^ with four 
insertions ip 2 (qi)/2, ip 2 (—qi)/2, ip 2 (q2)/2, <y9 2 (— g 2 )/2 one needs for calculation of the averaged 
asphericity ratios A^ AW , A^f 1 defined by (I4.21[) . The Green function G^' 4 \u) (see figure[5]) reads: 

G^\u) 



TR v ' (Q?+i)(g 2 2 +l) (^+i)(^I+l)[(gi+g 2 ) 2 +i] " (gf+l) 2 (g|+l) 2 
8/ u[2Ii(0,gi)+2Ji(0,<&)] u[2Ji(0 > ?i)+2Zi(0,g 2 )+27 2 (0 1 gi,g 2 )] 



3 l(gf+i)(gi+i)[(gi+g2) 2 +i] (g 2 +i)(g 2 2 +i) 

2uJi(0,gi) 2uJi(0, g2 ) l M/ 2 (0,0,g 2 ) l M/ 2 (0,0,gi) 

(g 2 +l) 2 [(gi+g 2 ) 2 +l] (g 2 2 +l) 2 [(gi+g 2 ) 2 +l] 2 (g 2 +l) 2 2 (g 2 2 +l) 2 
u[I 2 (0, 0, gi)+7 3 (0, 0, gi, g 2 )+7 3 (0, gi, g 2 , gi+g 2 )+7 3 (0, gi, gi, gi+g 2 )] 

(g|+l) 

m[/ 2 (0,0, g 2 )+7 3 (0,0, gi,g 2 )+/3(0,gi,g 2 ,gi+g 2 )+/3(0,g 2 ,g 2 ,gi+g 2 )] 

(gf+1) 

M[/ 2 (0,gi,gi+g 2 )+/ 2 (0,g 2 ,gi+g 2 )] _ M[/ 2 (0,gi,gi+g 2 )+/ 2 (0,g 2i gi+g2)] 

(g 2 2 +l)[(gi+g 2 ) 2 +l] (g 2 +l)[(g 1+ g 2 )2+l] 
u[27 4 (0, 0, 0, gi, g 2 )+27 4 (0, 0, gi, g 2 , gi+g 2 )+7 4 (0, 0, gi, gi, gi+g 2 )+7 4 (0, 0, g 2 , g 2 , gi+g 2 )] 



4u7i(0,0) 



(g 2 + l)(g 2 + l) (ql + l){ql + l)[(q 1+ q 2 y + l] [qf+lf 



(Bl) 



The Green function Gp (u, u>) reads: 
G^ 4) («,«,) = 



(g 2 +l)(g 2 2 +l) (g 2 +l)(g 2 2 +l)[(gi+g 2 ) 2 +l] (g 2 +l) 2 
«[2Ji(0, gi)+27i(0, g 2 )] u[2Ji(0, ?i)+2J 1 (0, g 2 )+27 2 (0, gi, g 2 )] 



3 l(g 2 +l)(g 2 2 +l)[(gi+g 2 ) 2 +l] (g 2 +l)(g 2 2 +l) 
w[2Ji(0,gi)+2Ji(0,g 2 )] w[2 Ji(0, gi)+2Ji(0, g 2 )+2J 2 (0, g u q 2 )\ 

' («? + l)(^ + l)[(«l+«2) a + l] (g 2 + l)(g 2 + l) 

2MJi(0,gi) 2uJi(0, g2 ) _ 1 ^7 2 (0,0,g 2 ) _ l M J 2 (0,0, gi ) 

'(g 2 +l) 2 [(gi+g 2 ) 2 +l] (g 2 2 +l) 2 [(gi+g 2 ) 2 +l] 2 (g 2 +l)2 2 (g 2 +l)2 

2wJi(0,gi) 2wJi(0, g2 ) l wJ 2 (0,0,g 2 ) 1 wJ 2 (0,0, gi ) 

"(g 2 +l) 2 [(gi+g 2 ) 2 +l] (g 2 2 +l) 2 [(gi+g 2 ) 2 +l] 2 {q 2 + l) 2 2 (g 2 2 +l)2 
m[7 2 (0, 0, gi)+7 3 (0, 0, q 1; g 2 )+7 3 (0, gi, g 2 , gi+g 2 )+7 3 (0, gi, gi, gi+g 2 )] 

(g 2 2 +l) 
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w[J 2 {0, 0, gi)+ J 3 (0, 0, gi, g 2 )+ J 3 (0, gi; g2, gi+g 2 )+ J 3 (0, gi, gi, gi+g2)] 

(gf+1) 

U[J 2 (0, 0, g 2 )+^3(0, 0, gli g2)+^3(0, gl, 92, gl+g2)+^3(0, 92, g2, gl+gg)] 

(gf + 1) 

w[J 2 (0, 0, g 2 )+ J 3 (0, 0, gi, g 2 )+ J 3 (0, gi, g 2 , gi+g 2 )+ J 3 (0, g 2 , g2, gi+g2)] 

(g 2 +l) 

u[/ 2 (0,gi,gi+g 2 )+/ 2 (0,g 2 ,gi+g 2 )] w[/ 2 (0,gi,gi+g 2 )+/ 2 (0, g 2 ,gi+g 2 )] 



(g 2 2 +i)[(gi+g 2 ) 2 +i] (g?+i)[(gi+g 2 ) 2 +i] 

w[J 2 (0,gi,gi+g 2 )+J 2 (0,g 2 ,gi+g 2 )] _ w[J 2 (0, gi, gi+g 2 )+ J 2 (0, g 2 , gi+g 2 )] 

(g 2 2 +i)[(gi+g 2 ) 2 +i] (g 2 +i)[(gi+g 2 ) 2 +i] 
m[2J 4 (0, 0, 0, gi, g 2 )+2/ 4 (0, 0, gi, g 2 , gi+g 2 )+/ 4 (0, 0, gi, gi, gi+g 2 )+^4(0, 0, g2, g2, gi+g 2 )] 

2 

w{2 J 4 (0, 0, 0, gi, g 2 )+2 J 4 (0, 0, gi, g 2 , gi+g 2 )+ J 4 (0, 0, gi, gi, gi+g 2 )+ J 4 (0, 0, g 2 , g 2 , gi+g 2 )] 1 



7 



4[mJi(Q,0) -iwJi(0,0)] 



(g 2 +l)(g 2 +l) (g 2 +l)(g 2 +l)[(g 1+ g 2 ) 2 +l] (g 2 + l) 2 (g 2 +l) 2 

(B2) 



The one-loop integrals Ii and Ji are explained in appendix A. 
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OopMM MaKpoMO/ieKy/i y xopouiux po3HHHHHKax: 
niflxifl TeopeTMKO-no/ibOBoV peHopMa/ii3aL4iMHoY rpynw 

B. BriaBaL^bKjR K. cdoh 0ep6erj2E] K). ro/ioBa^ 

I HCTMTyT 4>i3HKH KOHfleHCOBa Hwx cwcTeM HAH yKpa'mi/i, By.n. I. Cbeh L4i L4bKoro, 1, 79011 JlbBiB, YKpaiHa 
floc/iiflHMLibKnw L4eHTp npwK/iaflHoT MaTeMaTHKH, Vnii Bepci/ueT KoBeHTpi, CV1 5FB KoBeHTpi, Anrmn 
TeopeTMHHa cpi3HKa no/iiMepiB, VHiBepci/iTeT 0paR6ypry, D-79104 <J>paw6ypr, HiwieHHWHa 

y CTaTTi mm noKa3yeMO, sikhm huhom MO)KHa 3acTocyBaTH MeTOfl TeopeTHKO-nojibOBoT peHopMa/ii3aL4iMHoY 
rpynn fl/ia aHa/ii3y ymBepcajibHi/ix B/iacTHBOCTeCi cpopM flOBri/ix rHynKnx no/iiMepHwx /ia h L4K>riB y nopi/icTOMy 
cepeflOBninj. flo L4boro nacy TaKi aHajiiTHHHi po3paxyHKH b ocHOBHOMy TopKa/iwcb noKa3HM KiB CKewjiiHry, U40 
BM3HaHatoTb KOHCpopwia 14m Hi B/iacTM BOCTi no.niMepHi/ix MaKpoMO/ieKyji. ripoTe, icHyKDTb i/i iHLui cnocTepe>KyBa- 
Hi Be/iwHWHH, mo, 51 k i noKa3HM km CKen/iiHry, e yHiBepcajibHi/iMH (to6to He33/ie>KHM M H Bifl xiwiiHHoT CTpyKTypi/i 
hk MaKpoMO.neKy.fi, TaK i po3HHHHHKa), a OTwe MO>KyTb 6yTH npoaHajii30BaHi b Me>Kax niflxofly peHopMa/ii- 
3aL4i£iHoT rpynn. Mw L4i KaBnMocb nwTaHHfiM, akoT cpopMH Ha6yBae flOBra myHKa no/iiMepHa MaKpoMO.neKy.na 
y po3HHHi b npwcyTHOCTi nopi/icToro cepeflOBi/im,a. U,e nHTaHHa e cyTTeBUM p/\n po3y m i h H5i noBefliHKH Ma- 
KpoMOJieKy/i y KO/ioTflHux po3HMHax, n06.n1/13y MiKponopucTMX MeM6paH, a Ta ko>k y K/iiTM HHOMy cepeflOBumi. 
M1/1 po3nn:RflaeMO 3anponoHOBaHy paHiiue MOflejib no/iiMepa y cZ-BHMipHOMy npocTopi [V. Blavats'ka, C. von 
Ferber, Yu. Holovatch, Phys. Rev. E, 2001, 64, 041102] y cepeflOBumj is CTpy KTy pH n m m HeoflHopiflHOCTnMH, 
U40 xapaKTepn3yK)Tbca napHOK) Kope.n5mii7iHoio cpyHKujeio h(r), swa cnaflae i3 BiflCTa hhhd r 3riflHO CTeneHe- 
Boro 3aKOHy: h(r) ~ r~ a . 3acTocoByeMO niflxifl TeopeTMKO-no/ibOBoT peHopMa/ii3ai4i£nHoV rpynn i oujHioeMO 
BiflHoujeHHiR po3MipiB (-f?j?) / {Rq) Ta accpepw h HicTb A c i flo nepworo nopnflKy e = 4— d, 8 = 4— a-po3K/iafly. 

KyiiOHOBi c/iOBa: non'iMep, 3aMopo>KeniAiA 6e3/iafl, peHopMam3au,ii/iHa rpyna 
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